THE INTEGRAL/SPI RESPONSE AND THE CRAB OBSERVATIONS 
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ABSTRACT 



The Crab region was observed several times by IN- 
TEGRAL for calibration purposes. This paper aims 
at underlining the systematic interactions between 
(i) observations of this reference source, (ii) in-flight 
calibration of the instrumental response and (iii) the 
development and vali dation of the analysis tools of 
the SPI spectrometer (Ve drenne et al. 2003). It first 
describes the way the response is produced and how 
studies of the Crab spectrum lead to improvements 
and corrections in the initial response. Then, we 
present the tools which were developed to extract 
spectra from the SPI observation data and finally a 
Crab spectrum obtained with one of these methods, 
to show the agreement with previous experiments. 
We conclude with the work still ahead to understand 
residual uncertainties in the response. 

Key words: INTEGRAL/SPI; Crab; instrumental 
response; calibration. 



1. CRAB OBSERVATIONS 



During its first year in orbit, INTEGRAL observed 
the Crab region (the nebula and its pulsar) twice : 
in February (revolutions 39-45) and August (revo- 
lutions 102-103) 2003. The total observation time 
yields 1,723 ks (including annealing periods), per- 
formed in different modes : on- or off-axis, in staring 
or dithered mode, with a 5 x 5 or an hexagonal dither 
pattern. Because of its high luminosity and apparent 
lack of variability, this source is used as a reference 
for the spectrometer calibration. 



the redistribution matrix files (RMFs) are di- 
vided into three components describing (1) the 
events of the full-energy peak, (2) the Compton 
events interacting first in the detectors, (3) those 
interacting first in the passive material. They do 
not depend on the direction or the detector but 
only on the energy. 

the image response files (IRFs) describe, for a 
given input photon energy, the effective area of 
each detector for all directions in the field of 



Production The INTEGRAL/SPI instrument re- 
sponse is produced using a suite of Monte Carlo 
simulation software developed at NASA/GSFC 
based on the GEANT-3 package called MGEANT 
l|Sturner et al. 2003') . This production also required 
the development of a detailed computer mass model 
for SPI. 
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2. INSTRUMENTAL RESPONSE 



Response decomposition The response is de- 
composed into IRFs and RMFs : 



Figure 1. SPI photopeak effective area : comparison 
between ground calibration measurements (dots with 
error bars), and two successive IRF releases (dotted 
and plain lines). 



Absolute determination of the camera photo-peak ef- 
fective area was obtained from measurements (points 
in Fig.^ performed just before the launc h during an 
extensive ground calibration campaign IjAttie et al.l 
|2003; Schanne et al. 2003) of SPI in 2001 April at the 
Bruyeres-le-Chatel (BLC) site of CEA. These mea- 
surements were compared to simulations and led to 
the initial release of the SPI response (Fig. Q] dotted 
line) . 



Correction Initial analysis of ground calibration 
data was performed only on lines with energies at or 
above 60 keV. After the launch, it was found through 
analyses of early Crab observations, that the Crab 
flux below these 60 keV was underestimated com- 
pared to the spectrum expected. It enabled us to 
point out an over- estimation of the low energy effi- 
ciency. 

Additional analyses of 4 low-energy lines from the 
BLC data were performed. The 20.80 and 26.35 keV 
lines of the ^^^Am calibration source had initially 
been ignored because of their low statistics. The ^ 
30.8 and ~ 35.07 keV hues of ^^SBa, blended by the 
Compton component of higher energy lines, had also 
been rejected. A new response, corrected below 60 
keV by taking into account these 4 lines, was derived 
(Fig. ^ plain line) , without any assumption on the 
Crab spectrum. 



fhttp: //isdc .unige . ch/ index. cgi?Sof t+softj), 
adjusts the intensity of sources simultaneously 
with the scaling factors to apply to the back- 
ground model chosen. The best solution is 
searched for in each energy bin successively. 
The off-diagonal terms of the instrumental 
response are not taken into account in the 
deconvolution - all events are treated as pho- 
topeak events - and the resulting spectrum 
is in pseudo-photons. A ^^Spiros dedicated" 
redistribution matrix was derived using Monte 
Carlo simulations of both the SPI instrument 
and the Spiros software. This matrix can 
for example be used by the spectral-fitting 
program XSPEC, to fit a model to the real 
photon spectrum ; 

X SPEC 12. developped at NASA/GSFC 
l|Shradeill200^ is a new release of XSPEC in- 
cluding SPI specific packages - to be delivered 
soon. While XSPEC 11 can be used to fit a 
mathematical model to an already background 
substracted and deconvolved SPI source spec- 
trum, XSPEC 12 can work directly from the in- 
dividual detector spectra. Using, unlike Spiros, 
the full intrumental response ~ IRFs and RMFs 
- reconstructed for each source, it adjusts simul- 
taneously the model parameters of each of the 
sources considered and the scaling factors to ap- 
ply to the background model. 



3. DECONVOLUTION METHODS 




Figure 2. Count rates in the 19 Ge detectors in the 
20-8000 keV hand, during a Crab on axis pointing. 
From left to right : total count rate, estimated hack- 
ground count rate and their difference. The pattern 
of the projection of an on axis source through the 
mask on the detector plane, with its 120 degrees sym- 
metry, is visihle. 

Several tools have been developed to extract source 
spectra from SPI observation data, including Spiros 
and XSPEC 12. 



in its spectral mode, Spiros IjSkinner fc ConnelJ 
l200l . available in the INTEGRAL off- 
line scientific analysis (OSA) distribution 



To account for initial discrepancies in the Crab 
derived flux between Spiros and XSPEC 12, a 
new version of the ^^Spiros dedicated" redistribu- 
tion matrix mentioned above was recently derived 

at NASA/GSFC. 

However, the main challenge is to properly model the 
background which represents a huge fraction of the 
total measured count rate (Fig. O and 3). 
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Figure 3. Total count rate (upper line) and estimated 
background count rate (lower line) during a Crab ob- 
servation : the signal of interest represents only a 
tiny fraction of the total count rate. 
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Figure 4- SPI Crab spectrum extracted from 44 NRT science windows in INTEGRAL revolution 44^ using OS A 
3.0 and a saturated Ge background model fitted detector by detector with Gaussian statistics. 



Spectral studies of Crab and other reference sources 
were done with both Spiros and XSPEG 12, enabhng 
us to cross-cahbrate these tools. The resulting spec- 
tra are now consistent in shape and normalization. 
An initial discrepancy in normalization was reduced 
after a new redistribution matrix was derived for 
Spiros. Some differences remain which are still under 
study. 



4. CRAB SPECTRUM 

The Crab spectrum in Fig. 0] was extracted using as 
standard a method as possible, so that any scientist 
undertaking to analyze SPI observations shall be able 
to reproduce such a spectrum. It corresponds to the 
first 44 near real time science windows of revolution 
44, performed with a 5 x 5 dither pattern. Distribu- 
tion 3.0 of the observation analysis software (OSA) 
was used and the spisciencc-analysis script was run 
mainly with default parameters, the most important 
of which concern the background : spiback produced 
a background model based on the count rates in sat- 
urating Ge detectors and Spiros adjusted the source 
flux and background count rates in each energy bin, 
assuming the time variation given by the model and 
computing the best detector to detector ratios. 



Table 1. Flux (ph/crr? /s) and power law index of 
total Grab gamma-ray emission. 



Instrument 


Index 


50-100 keV flux 


OSO-8 

GRIS 

CGRO/OSSEt 

CGRO/BATSEt 

SAX/PDSt 

INTEGRAL/SPI 


2.00 ±0.06 
2.15 ±0.03 

2.19 ±0.03 

2.20 ±0.01 
2.13 ±0.01 
2.17 ±0.01* 


6.4110-2 
4.5210-2 
5.6810-2 
6.8310-2 
4.9210-2 
7.08 ±0.03* 10-2 



t private communication 
* statistical error only 



Above 1 MeV, the flux uncertainty grows up as the 
source flux gets smaller and the basic background 
model used starts to show its limits : the features vis- 
ible above 1 MeV correspond to instrumental back- 
ground residuals. 



A power law F{E) ~ K (j-^) with a spectral 
index a = 2.169 ± 0.008 and a normalization K = 
14.44 ± 0.44 ph/cm^/s/keV fits well to the spectrum 
(x^/d.o.f = 1.31) in the 40 keV to 8 MeV energy 
range. 

Imposing the canonical photon spectral index of 
2.10 found by X-ray experiments for the entire 
Crab, we estimate (x^/d.o.f = 2.39) a normaliza- 
tion of 10.80 ±0.03 ph/cmVs/keV at 1 keV, to com- 
pare with a value of 9.59 ph/cm^/s/keV given by 
iWillins-ale elTTI lf200l . 

A broken power law fits only slightly better (x^/d.o.f 
= 1.29) than a single power law. Imposing a low 
energy index ai = 2.10, we find a high energy index 
a2 = 2.19 ± 0.01 with a break around 61 ± 6 keV 
and a normalization of 11.03 ± 0.05 ph/cm^/s/keV 
at 1 keV. 

Using a different energy range or background han- 
dling method, the fit parameters found vary slightly. 
Even without taking into account systematic uncer- 
tainties due to the calibration and background han- 
dling works still in progress, Table Q] shows a very 
good agreement between SPI and other gamma-ray 
experiments. 
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5. CONCLUSION 



After a few months of fine tuning of both the in- 
strumental response and the deconvolution tools, the 
Crab spectra extracted from INTEGRAL/SPI obser- 
vations with the various software available are com- 
patible with each other and rather consistent with 
previous experiments. 

Some further work will be necessary to build more 
elaborate background models and become more con- 
fident in the fluxes derived. Although the response 
and the extraction software might still evolve in the 
future especially to account for the loss of detector 
2, the tools available today are sufficient to study 
point sources observed with INTEGRAL. 
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